1 



Resonances and the thermonuclear reaction rate* 
O ■ M. S. Hussein'^, M. Ueda^, A. J. Sargeant'' and M. P. Pato'^ 

o ■ 

■ ^Nuclear Theory and Elementary Particle Phenomenology Group, Instituto de Fisica, 
Universidade de Sao Paulo, Caixa Postal 66318, 05315-970 Sao Paulo, SP, Brazil 



(N 



^Akita National College of Technology, lijima Bunkyo-cho 1-1, Akita, 011-8511, Japan 



>■ ' 1. Abstract 

cn ■ 

^ ' We present an approximate analytic expression for thermonuclear reaction rate of 

■ charged particles when the cross section contains a single narrow or wide resonance de- 
O . scribed by a Breit-Wigner shape. The resulting expression is uniformly valid as the 

! effective energy and resonance energy coalesce. We use our expressions to calculate the 
reaction rate for ^^C(p,7)^^N. 

O I 2. Introduction 

Evaluations of thermonuclear reaction rates require the folding of nuclear cross sections 
with the Maxwell-Boltzmann distribution. In recent compilations of such rates for charged 
^ ' particles [Jj (^=1-14), |2j (A=20-40) a combination of numerical integration and analytic 

■ techniques appropriate to the energy dependence of the cross section for each reaction 
pair were used. A review of the analytical techniques of nuclear reaction rate theory 
can be found in Ref. [3]. Ref. ^ has obtained results different to those obtained in 
Ref. due to the inclusion of resonances in their numerical integrations. Ref. PP 
calculated the contribution of narrow resonances in the simplest possible manner which is 
to approximate the Maxwell-Boltzmann distribution by its value at the resonance energy. 
This approximation can be expected to be good for resonances which do not overlap 
significantly with the Gamow peak. 

In Ref. jS] we developed an asymptotic expansion for the thermonuclear reaction rate 
in terms of the effective astrophysical ^-factor, Ses, using the method of Dingle [Hj. The 
method may be used in cases where S{E) can be reliably expanded as a Taylor series. 
Two alternative expressions Scs where obtained by expanding S{E) about E=0 and about 
E=Eq, Eq being the effective energy of the Gamow window. From these expressions, all 
the approximate formulae for Ses commonly used [7] when S{E) is slowly varying may be 
obtained as special cases. The validity of the expressions derived in Ref. j5j is limited by 
the radii of convergence of the Taylor series expansions of S{E), that is, by the location 
of the poles of S{E). The poles may be due bound states of the composite nucleus. 



* Supported by FAPESP. 



2 



Illustrative is the case of ''Be + p — + 7 jHI for which S{E) has a pole at E=—Eb, 
£■5=137.5 keV being the binding energy of ^B. (Ref. discusses effects other than a 
sub-threshold pole which can produce the low energy rise seen in the S'-factors of several 
capture reactions.) There may also be poles in S{E) due to resonances of the reaction 
pair in which case the cross section may be parameterised in the region of the resonance 
by the Breit-Wigner form (see Eq. (j21) below) which has two poles, or the Lorentzian form 
1101 which 
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Figure 1. The reaction rate for ^^C(p,7)^^N as function of Tg. The circles show the rates 
from Ref. P, the squares Ref. |TI], the solid lines the uniform approximation, Eq. (jSJ, and 
the dotted lines the resonance-plus-tail approximation, Eq. (jl)). We took the resonance 
parameters from Ref. ^| (£^=457 keV and ^ Fp(i?r)=39 keV [both being laboratory 
frame values] and a (Er) =130 ji b). 



For a given reaction cross section (t{E), the Maxwellian averaged reaction rate per 
particle pair is given by [TB] 
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where is the reduced mass of the two collision partners, ks is the Boltzmann constant, 
and T the environmental temperature. Here, we consider a (p, 7) resonant reaction for 
which the cross section is given by 



a(£) = ——LU 



T,iE)T, 



2/i £ {E- Er^ + VII ^ 



(2) 
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where a; is the statistical spin factor, is the 7-decay width and Tr is the total width 
at the resonance energy. The energy dependence of the proton-decay width Tp{E) is 
determined by the Coulomb barrier for low enough energies and may be written as 
Tp{E)='ypexp [—27Tri{E)], with 7p energy independent. The Sommerfeld parameter, ri{E), 



is given by 2Tirj{E) = ^ Eg / E , where EG=b'^=2fic^ {■nZiZ2(y) , a being the fine structure 
constant. The effective energy of the Gamow window at temperature T is found to be 

^0 
Su 



E]!\bT/2\^'\ 

Dstituting Eq. Q into Eq. and introducing the dimensionless variables x=E /ksT , 



Xr=Er/kBT and xi=Tr/2/kBT one can rewrite the reaction rate as 

' TT/i j KKbT J 2/i Xi Jo X — Xp 



where F{x)=x+yxG/x, Xp=Xr+ixj and xg=0'^=Eg/ ksT . The function F{x) has a sin- 
gle minimum in the the interval [0, 00) from which it increases monotonically to in- 
finity as X approaches zero and infinity. The position of the minimum is given by 

r 1/2 1 ^/'^ 

XQ=EQjkBT= Xq /2 . The principal contributions come from the neighbourhood of 
the critical points xq and Xp j6fl4j . 

3. Conventional approximation 

When the resonance and the Gamow window are well separated the contribution from 
each to the reaction rate can be estimated separately and the two contributions summed. 
Let us write {av) ~ {(tv)^ + {(tv)gi where {(Tv)r is the contribution from the resonance 
and (cTf )(7 the contribution from the Gamow window. The resonance contribution may be 
estimated j[13i by approximating the exponential in Eq. Q by it's value at the resonance 
energy, Xr, so that dxe~^'^''y[{x - Xrf + x]] ^ g-^v-^^^cAv (^^^2 + Or) jxi, where 
0r=tan~-^ [xrjxi). 

Similarly, the Gamow window contribution is estimated by approximating the Breit- 
Wigner factor by it's value at the location of the Gamow peak, xq, so that /q°° dxe'^'^^^ I\{x— 
Xr)^ + Xj] ~ |v^7rre^'^/[(xo — x^)^ + x|]. We have also used the Gaussian approxima- 
tion F(x)-F(xo)^(x - xo)^ / HHj where A^=y/8/F"(xo)=4^xo/3 is the effective 
width of the Gamow window in units of ksT and r=F(xo)=3xo. 

On thus obtains 



{ov) ^ (2//i)^/^— ^ ^7rV2^^(£;,)e-^'■/^.T-y^^^^(^^)^-. 



(4) 



where we have introduced the astrophysical S'-factor, S{E)=a{E)Eev^^^ and the width 
of the Gamow window A=4(£'o/csT/3)^/^. The second term of Eq. (0} is the zeroth order 
term of an asymptotic expansion of Ig in powers of 1/r jS] 



4. Uniform approximation 

We introduce a new integration variable t implicitly through F{x)—F{xo)=t'^ and 

dt X- 



rewrite the reaction rate in terms of the function $(t)=%-^z^- The position of the 
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pole in t-space is then given by tp 



1/2 



In order to obtain an expan- 



'Xc/Xp - T 

sion which is uniform made to approach or coincide with xq, $(t) is expanded as 

follows mHaiE]: $(t)=E^=o [«n + it- tp)/3n] [t{t - tp)]". This expansion has the prop- 
erty that it gives the value of the function and all it's derivatives exactly at both 
t = and t = tp. Let us approximate by the truncation $o(i)=«o+(^ ~ tp)Po- The 
expansion coefficients are found to be ao=l and Po=l/tp+2y/T/3/ (xq — Xp). Introducing 
the complementary error function |T7] we obtain 



{av) ^ (2//i) 



1/2 



(fcBT)3/2 



TS{Er)'^ \ i^e-'hrici-itp) + - + AS{Eo 



(5) 



5. Reaction rate for ^^C(p,7)^^N 



For the 5'-factor at the resonance energy we use S{Er)=^^^^^e''/^^^ and at the 
effective energy S{Eo)=S{Er) (^Eo-e )^+r^/4. ' "where the resonance strength is given by 
(a;7)r=^^^^^^^. We have calculated the reaction rate for ^^C(p,7)"'^^N using Eqs. (0) 
and (jH) including only the lowest resonance and the results are displayed in Figure ^ The 
exact integral is not shown as it is indistinguishable from the uniform approximation for 
this case. It is apparent that the resonance-plus-tail approximation cannot be used at 
temperatures for which E^r^E^. For comparison we have also shown the rates of Ref. [TT] 
which used a different analytic approach and those of Ref. which also included the 
contribution of the second resonance at 1698 keV. 
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